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Oh! Abstract 

We consider the two-Higgs-doublet model with exphcit CP vio- 
lation, where the effective Higgs potential is not CP invariant 
at the tree-level. Three neutral Higgs bosons of the model are 
5^ \ the mixtures of CP-even and CP-odd bosons which exist in the 

CP conserving limit of the theory. The mass spectrum and tree- 
level couplings of the neutral Higgs bosons to gauge bosons and 
fermions are significantly dependent on the parameters of the 
Higgs boson mixing matrix. We calculate the Higgs-gauge bo- 
son, Higgs-fermion, triple and quartic Higgs self-interactions in 
the MSSM with explicit CP violation in the Higgs sector and 
CP violating Yukawa interactions of the third generation scalar 
quarks. In some regions of the MSSM parameter space substan- 
tial changes of the self-interaction vertices take place, leading to 
significant suppression or enhancement of the multiple Higgs bo- 
son production cross sections. 
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1 Introduction 



General interest to the models with two (and more) Higgs doublets is main- 
tained by the absence of a convincing argument in favor of only one genera- 
tion of Higgs bosons when there are three known generations of fundamental 
fermions. Models with extended Higgs sector provide richer physical possi- 
bilities than the standard scheme with one doublet. One of them is the pos- 
sibility to introduce CP violation beyond the Cabibbo-Kobayashi-Maskawa 
(CKM) mechanism, by means of the Higgs boson exchange amplitudes with 
complex Higgs boson-fermion vertices. Complex couplings can be generated 
either spontaneously [1], when the vacuum expectation values of the Higgs 
fields are complex and couplings of the CP invariant tree-level Higgs poten- 
tial are real, or explicitly inserted [2] on the level of SU{2) x t/(l)-invariant 
potential terms, when the complex vacuum expectation values of scalar fields 
correspond to the minimum of hermitian potential with complex couplings, 
which is not CP-invariant (CP-invariance softly broken by the mass terms). 

Various representations of the SU{2) x [/(l)-invariant two-doublet Higgs 
potentials have been considered in the literature. The two-doublet models 
with spontaneous CP violation [1, 3] make use of the potential of general 
structure —fj.'^f^ + A</7^ without the dimension two /i^2"terms. Models with 
explicit CP violation use either the potential with trivial minimization ([2], 
see also [4]) or the potential with complex coupling /i^2 of the dimension 
two terms and complex couplings A5, \q and A7 in front of the dimension 
four potential terms [5, 6], similar to effective potential of the minimal su- 
pesymmetry (MSSM). Standard transformation (diagonalization) procedure 
from the level of primary fields which are the components of scalar doublets 
in SU{2) X t/(l)-invariant potential terms {SU{2) x U{1) eigenstates), to 
the physical fields (mass eigenstates) of Higgs bosons should be consistently 
performed to respect the SU{2) x C/(l) invariance and the minimization of 
the potential. Wc consider the diagonalization for the two different two- 
Higgs-doublct potential forms in much details. A special case of the general 
two-Higgs-doublet model is represented by the Higgs sector of MSSM. Sub- 
stantial radiative corrections to the Higgs boson masses and couplings are 
induced at the mz scale mainly by the third generation quarks h and the 
third generation scalar quarks [7]. In the special case of MSSM the multi- 
parameter space of general two-Higgs-doublet model is significantly reduced, 
providing possibilities of much less ambiguous phenomcnological predictions. 

Phenomenological consequences of the CP violating Higgs-third genera- 
tion squark Yukawa interactions in the Higgs-fermion and the Higgs-gauge 
boson sectors have been considered in [5]. Wc focus mainly on the self- 
interactions of Higgs bosons. Experimental observation of the scalar boson 
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signals should be followed by the verification of Higgs mechanism as the 
essence of the gauge boson and fermion mass generation. Self-interactions of 
the Higgs fields lead to untrivial structure of the vacuum state with nonzero 
(and possibly complex) field tensions, initializing the spontaneous breakdown 
of SU{2) X U{1) symmetry. Reconstruction of the Higgs self- interaction 
potential from the data on multiple (mainly double and triple) Higgs bo- 
son production cross sections [8] requires the experimental measurements of 
triple and quartic Higgs boson self-interaction vertices, which is nontrivial 
but valuable task for a future high luminosity colliders, such as LHC and 
TESLA. 

In section 2 we discuss the diagonalization of CP invariant Higgs potential, 
represented in two different forms, in the general two-Higgs-doublet model 
(THDM) and consider the special case of the Higgs sector in the minimal 
supersymmetry (MSSM). In section 3 we introduce complex parameters of 
the SU (2) X U (1) invariant potential terms and discuss the diagonalization of 
the THDM potential which acquires explicit CP violation. In sections 4 and 
5 we calculate the Higgs-gauge boson, Higgs-fermion and Higgs self- couplings 
in the MSSM with CP violation. 



2 Diagonalization of the mass matrix in the 
general two-Higgs-doublet model 

Two representations have been used for the two-doublet Higgs potential. The 
first representation [2, 4] 

V{^U ^2) = XMVl - if + \2i^t^2 - if (1) 

+A3[(v^^V'i - f ) + ivtv2 - i)? 

+ M{vtVl)i'^2V2) - (V^^V^2)(<^2<^l)] 

+X,[Re{ipt<P2) - ^Re(e*«)]2 + Xe[lm{iptip2) - ^Im(e*€)]2 
where Aj are real constants and the SU (2) doublets (pi^2 have the components 

(pi = {-iwf,-^{vi + hi + izi)}, (p2 = {-iw2,-^{v2 + h2 + iz2)}. (2) 

wi^2 are complex fields and 2:1 2, /ii,2 are real scalar fields. At positive Ai, ...Ae 
each term of the potential V{Lpi,Lp2) is obviously positive and its zero min- 
imum is achieved if the vacuum expectation values of < (^i >, < (^2 > are 
taken in the form 

< fi >= -^{0, vi}, < f2 >= -^{0, V2e'^} (3) 
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In the case of A5 = Ag (corresponding to the CP conserving MSSM-hke 
potential, see below) the last two terms in (1) form the modulo squared, 
and the phase ^ can be removed from the potential by the U{1) rotation 
of (f2, that does not affect minimization. In this section we will consider 
the case A5 7^ Ag, C =0- Substitution of (2) to (1) gives a bilinear form of 
the mass term with mixed components Wi, hi, Zi, which can be diagonalized 
by an orthogonal transformation of the fields in order to define the tree 
level masses of Higgs bosons. In the CP conserving case the potential terms 
involving zi, fields from the imaginary parts of </7i, (/?2 doublets and h\, 
fields from their real parts do not mix, so the mass terms are diagonalized by 
a separate two-dimensional rotations of the z\, Z2 and the hi, h2 fields. The 
resulting spectrum of scalars consists of two charged H"^, three neutral h, H, 
scalar fields, and three Goldstone bosons G. This procedure is described 
in many papers (for instance, [4, 9]). The Wi^2 sector is diagonalized by the 
rotation of Wi,W2 ^ H,G 

wf = -H^S0 + G^cp, wf = H^cp + G^sp (4) 

defined by the angle 

tg/3 = ^ (5) 

and leading to the massless G field and the field of massive charged Higgs 
boson H^, nT'H± — ^4(^1 +f2)/2- The Zi^2 sector is diagonafized by the 

rotation Zi, Z2 ^ A'^jG defined by the angle f3 and giving again one massless 
field G and the field of CP-odd Higgs boson A^ with the mass m\ = \Q{vf + 
w|)/2. Finally, the hi, /i2 sector is diagonalized by the rotation hi, h2 ^ h,H 
defined by the angle a 

. r, 2mi2 „ mil - m22 

^ {mil - m22y + 4mi2 y (mn - 77122)^ + 4mi2 

where 

mn = \[Avl{\i + \z) + vlW 
m22= i[4vi(A2 + A3)+w?A5] 
mi2^ |(4A3 A5)viV2 

giving two massive fields of CP-even Higgs bosons H, h with masses 



'^Hh^ "^11 + "^22 ± V (^11 - ^22)^ + 4m?2 (7) 
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The diagonal mass term of scalar fields and their triple and quartic self- 
interaction vertices can be explicitly obtained by the substitution of the fol- 
lowing expressions for Aj to the potential V{ipi,ip2) (1)^ 

2^^[5c-/3m^- J.._,m|,] + ||A5 (8) 

^3= 2^[-f?7^i + f?7<]-iA5 



2^2 1 S2l3"'h^ S2^"'Hi 4' 

6 = ^"^AO 



where we used the notation = vf + w|, Sq, = sina, Ca = cosa. Diag- 
onalization of the mass term takes place for arbitrary A5, which is a free 
parameter of the model. 

The second representation of the Higgs potential 

U{ipi,ip2)^ - Hl{(pt 'Pi) - f^lif 2^2) - f^i2{VlV2 + Vtfl) (9) 

+\l{iptipiy + \2{<P2^2f + >^2.{^t^l){^2^2) 
+X^{iptV2){vtV\) + T\(vtV2){vtV2) + ivtVllivtVl)] 

originates from the general SUSY action after the integration over Grassman 
variables and introduction of the soft SUSY- breaking terms (see [4]). It is 
easy to check that the potentials (1) and (9) are equivalent if parameters Aj, 
//f , ^2-1 1A2 Aj are related by the formulae 



Ai = Ai + A3, A2 = A2 + A3, A3 = 2A3 + A4, (10) 
^.Ag V _ A5 Ag 
2 ^ 2 ' ' ~ 2 2 



T _ . A5 Ag T _ A5 Ag 2 _ \ ^l'"2 

M — —M + Tr~'~~o~' '^s — t; TT, 1^12 — -^5 o 



and 



lij = Xivl + Xsvl + Xsvl, 111 = A2V2 + A3'yi + Xsvj (11) 

Unlike the potential (1) where the minimization is obvious, the symbolic 
structure of (9) does not demonstrate evidently its minimum. The substi- 
tution of (2) to (9) gives linear terms in the component fields zi^2, /ii,2 (or 
physical fields h, H, A) and unless some additional conditions to remove the 
linear terms are imposed, we are not in the minimum of the potential. So 
the equations (11) which set to zero the terms which are linear in compo- 
nent fields are the minimization conditions. The diagonalization of U {(fi, (^2) 
takes place for arbitrary /x^2 parameter. 
Inverse transformation (10) has the form 

\ _\ A4 A5 A5 X _T A3 A4 A5 A5 

M-X^------ + -. A2-A2-y-y-y + y (12) 
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As — ^ + ^ + ^ — A4 — — A4 — A5 + A5, Xq — — 2 A5 + A5 

so masses of the CP-even scalars and their mixing angle a (6), (7) in the case 
of potential U{ipi, 1P2) can be easily obtained using 



frill + rn22 = v^Xi + vl\2 H — —, mu - 777,22 = v^Xi - vl\2 - ctg 2/3/7-12) (13) 

82/3 

2/7712 = t'l^'2(A3 + A4 + A5) - /7i2- 

The diagonal form of U{ipi, 992) and the physical scalar boson self-interaction 
vertices are obtained by substitution of the following expressions for Aj and 
lii to (9): 

2j^[($)X + (5)'<-f/^?2] (14) 

M%fml + {f^rm%-^pil2] (15) 

A3= ^[2777^,-^ + ^(^771,- 777^] (16) 

^4= ^(1^+ 777^- 27771,,) (17) 

A5= Mj^-O (18) 

1^1^ Ai7;2 + (A3 + A4 + A5)f -/7?2tg/3 (19) 

lil^ A27;i + (A3 + A4 + A5)f -/7?2Ctg/3 (20) 

Conditions (14)- (18) ensure the diagonal form of the mass term expressed in 
physical fields h, H, A, and (19)- (20) are the minimization conditions. 
Two parametrizations for the Higgs boson self-interaction vertices can be 
used in THDM. In the first parametrisation [10] /7f2 is a free parameter and 
As is defined by (18). In the second one A5 is a free parameter and fil2 is 
equal to S/3C/3(w^As + rn^). Complete sets of Feynman rules (unitary gauge) 
for the triple (/7I2 and As parametrisations) and quartic {11I2 parametrisation) 
Higgs boson interactions in the general two-Higgs-doublet model with CP- 
conservation are shown in Tables 1-2 ^. In the case of MSSM potential at 
the scale Msusy (see (29)) As =0 and it follows from (8), (10), (11) that /7I2 
is fixed and equal to rn^spcp. 

Two additional terms of the dimension 4 can be constructed using the 
complete set of SU{2) x U{1) invariants (fitfi, <f2f2, Re(fi(f2 and lm(fi(f2 

^ These sets were obtained by means of the LanHEP package [11], see 
http://theory.siiip.msu.ru/~semenov/lanhep.html Six misprints in sign of the sec- 
ond term in the factor — s^sp) that have taken place in [10], are corrected in Table 3 
in the expressions for vertices hhhh, AAAA, hhAA, H+H-H+H', H+H~AA, H+R-hh. 
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(a detailed discussion of all possible potental forms can be found in [12]). 
These terms are usually added to the U{ipi, (^2) with the couplings Xq and A7 

U{if,,<f2) = U{<f,,<f2) (21) 

The diagonal form of U{(pi, ^2) at the local minimum takes place at arbitrary 
/x^2) Ae, A7 and can be achieved by means of the substitution (14)-(20) with 
additional Ae, Ay terms in the right-hand side: 



Ai = ^[(^)X + (^)^m|-^/.L] + 7(W/?-3A6tg/3) (22) 
A2 = iLiiv^'^l + &"^H - + ^(A6Ctg3/3 - 3Arctg/3) (23) 

A3 = i2m^.-^ + i^(m|-mD]-^ctg/3-^tg/? (24) 



1 / /^i2 , 2 o 2 \ Ae „ A 



A4 = -^{^ + m\-2mU)-'-^otg(3-'-^tg(5 (25) 



1 / A*i2 2 \ Ae ^ A 



A5 = - mi) - - (26) 

spcp 2 2 

2 22 

fil = Ai^;? + (A3 + A4 + A5)Y-M?2-tg/3+^(3A6Ctg/3 + A7tg/3)(27) 

2 22 
l4 = A2^;2' + (A3 + A4 + As)^ - //^ctg/? + ^(Agctg/? + 3A7tg/3)(28) 

Our expressions for the redefined A4 and A5 are the same as given in [13]. 

The potentials (1) and (9) can be reduced to the MSSM potential in some 
regions of the parameter space which we are going to discuss. The potential 
V{(pi,ip2) (1) has eight parameters: two vev's Vi, V2 and six couplings Aj 
(i=l,...6). Eight parameters of the potential f/(v'i,v?2) (9) fJ^i, /i2, ^^12 and 
Aj (i=l,...5) can be found using (10), (11). From the other side, in the Higgs 
sector we have eight physical parameters: the mixing angle /5 and ly-boson 
mass niw, mixing angle a, the parameter //f2 four masses of scalars m^, 
iTiA, mH±- The niw is fixed experimentally maintaining the constraint 
on the f 1, f2, v"^ = v\ + V2 = 4m^//e^ ■ slv?9y/ which follows from the Higgs 
kinetic term D^j^ipD^ip {g = e/sinO]^, 9y/ is the Weinberg angle). So the 
Higgs sector of THDM with the potentials (1) or (9) is described by a seven- 
dimensional parameter space. In the case of superpotential five additional 
constraints are imposed, relating all Higgs boson self-couplings Aj, (i=l,...5) 
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to the gauge coupling constants at the energy scale Msusy [14]: 

\SUSY _ \SUSY _ 9^ + 9l \SUSY _ 9^ ~ Ql \SUSY _ 9^ \SUSY _ n 
'^l - ^2 — g ) ^^3 — ^ ' '^i — 2"' 5 — 

(29) 

It follows that the four Higgs boson masses and the two mixing angles are 
defined by two independent parameters. One can choose, for example, the 

ri,r2 parametrization [15] (ri^2 = i^\h/'^%) well-known m^, tgf3 

parametrization. In order to reduce the general two-Higgs-doublct model 
vertices to MSSM at the scale Msusy it is convenient to use the a, {3 
parametrization: 

2 2 '^a+/3 2 2 2 2 ^2{a+l3) 2 2 

= m^C2/3 , rriH = m2C2/3 , = mz , A*i2 = m^s^cp. 

(30) 

Substitution of these expressions to the vertex factors in Tables 1,2 after 
trivial trigonometric transformations reduces them [10] to a simpler MSSM 
factors (see [4]). However, (30) are no longer valid at the energy scale mw 
where \f^^^ couphngs and masses of Higgs bosons are significantly changed 
by radiative corrections and the effective two-Higgs-doublet potential should 
be described in the complete seven dimensional parameter space. Practi- 
cal calculations of the radiatively corrected masses and/or couplings can 
be conveniently carried out using results of the two approaches, renormal- 
ization group (the HMSUSY package [16] or the analytical representation 
[17]) and diagrammatic (the FeynHiggsFast package, see [18]). Two different 
parametrizations can be used for these approaches. 

In the RG approach it seems convenient to use the two-Higgs-doublet 
model parameter space m^, tg/3, Ai, ...A5. In the following we shall take into 
account the Ae and Ay terms defined by (21), so the parameter space will be 
nine- dimensional. RG evolution of the coupling constants Aj from the energy 
scale Msusy to the electroweak scale mw defines the Ai,...A5 in (22)-(26) and 
the parameters Ag, A7. At a given m^, tg/9, Ae, Ay the parameter /i^2 ^-iid 
mH± are fixed by the conditions (25) and (26), the parameters ^\ and 
are fixed by (27) and (28), a, ruh and mn can be found using the equations 
(22)-(24). If we denote the deviation from the coupling Af^^^ at the MSSM 
scale by AA^ 

[SUSY \ \ A\ \SUSY 
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2(^1,2 ~ ''^1,2) — ^^1,2, A3 4 — A3,4 — AA3,4, — A5,6,y — AAs^ej 

we find the mixing angle (introducing the notation g\ -\- — g^rrr^/m 
9^-9\^ ^'(2 - m|/m^) while using (22)-(26)) 

^ g2/3(m^ + m|) + t;'((AA3 +_AA4)52/3 + 24AA6 + 24AAy) 

"02/3(771^ - m|) + t;2(AAi4 - AA2S^ - AA5C2^ + (AAg - AAy)s2/3) 
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CP-even Higgs boson masses and the parameter 



2 2 2,2 2 ron\ 

-v'(AAi4cJ + A\2slsl + 2(AA3 + AX^)caCf3SaSf} + AAsl^sJ + s^cj)) 

+2Sa+/3(AA6CaC/3 + AAyS^S^) 
2 2 2,2 2 

-t;^(AAis^cJ + AA24s5 - 2(AA3 + AX4)caC0Sc,sp + AAgls^sJ + 4cJ)) 

-2Ca+/3(AA6SaC/3 - AXjCaSp) 



2 



"i? — — 

mjj± = + - y (^'^s ~ ^-^4) (33) 

- - 

/i?2 = Sf3Cp[ml-— {2 AX5 + AXectgp + A\^tgf3)] 
with the minimization conditions 

IJ'l = ^"^|c2/3 - /x?2"tg/5 (34) 

t|2 _ _ _ _ _ _ 

- — [AAi4 + (^^^3 + AA4 + AAs)^ + 3AA6S^C;3 + AX^^] 
2 

9 3 

— — — — — C/q — 

- — [AA24 + (^-^3 + AA4 + AA5)4 + A^6— + 3AA7S/3C^] 
2 "^5/3 

These expressions can be straighforwardly used to calculate the radiatively 
corrected masses of Higgs bosons and the mixing angle a in the MSSM with 
the help of a solution of the RG equations for Ai, ...Ay. Apparently, in the 
RG approach Feynman rules in terms of Aj couplings are more convenient 
than rules in terms of Higgs particle masses. 

In the diagrammatic approaches to calculation of the radiatively corrected 
masses [18] the corrections to m^, in a and m^f± are extracted from the 
renormalized Higgs boson self-energies (usually radiative corrections to only 
the GP-even Higgs boson masses are calculated). The set of 7+2 independent 
parameters inherent for the diagrammatic approaches could be niA- tgf3, a, 
I1I2, iTih-i ^H, ^H±-, and Ae, A7. At a given m^, tg/5, Ae, A7 the /i^2 parameter 
can be fixed at the value rri^sgcp, and a can be calculated using the renor- 
malized self-energies correction [18] to the relation valid at the Msusy scale 
m\ + m\ = —S2a/s2f3{TrL\j — m'j^. Then A4 is defined by (25) and Ai,...A3 can 
be found using (22)- (24). In the diagrammatic calculations Feynman rules 
in terms of the radiatively corrected Higgs boson masses look more natural. 
Substitution of the radiatively corrected Higgs masses to the tree-level Higgs 
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vertex factors is expected to give results very close to those obtained from the 
loop corrections to Higgs vertex at the SUSY scale (see the discussion in the 
last ref. [8]). It has been shown in [19] on the example of hhh and hhhh ver- 
tices (and for the case of diagonal third generation squark mass matrix) that 
large radiative corrections to the vertex factors calculated diagrammatically 
can be absorbed in the radiatively corrected Higgs boson masses. 

Other parametrizations in the two-Higgs-doublet model are of course pos- 
sible, but they should be carefully introduced to respect the minimization 
and diagonalization conditions (22)- (28). The introduction of scalar particle 
masses and mixing angles inconsistent with them violates either diagonal- 
izaton of the potential or its SU{2) invariance, even if the minimization 
conditions remain valid. 



3 CP violation in the two-Higgs-doublet model 

CP transformation of the scalar doublet CPtpP^C^ =Ccp 'P~^ (the phase 
factor ICcpI =1) changes the sign of imaginary part lm{ipfip2) in the Xq term 
of potential (1), so if A5 7^ Xq, Xq 7^0 and ^ 7^0, CP symmetry is broken there 
explicitly. In other words, the dimension two terms of (1) appear with the 
complex parameter //fg 

^[(pt^P2 + vt'fi - ^^l^^2Cos{]^ + ^[-ii^tV2 - ^t^i) - Vi^asin^]^ (35) 

^ - ^)[{¥^tV2f + ivtvi?] + (y + y ) vt^2 <ptvi 
(A5COS4 — ^A6Sln4J (A5COS4 + ^A6Sln4j <^2 V'l 



so we find 

iA.2 = ^y^l'^scos^ - iAesin^). (36) 

CP is (softly) broken by the ^\2ViV2+ f^u vt^i terms. In this special case 
when the same phase ^ is involved both in the potential and the vacuum ex- 
pectation value of (/?2, diagonalization and minimization of CP transformed 
potential become less transparent. It is more convenient to analyse the po- 
tential form (21) for the general case of complex parameters with arbitrary 
phases. In the following we shall consider the hermitian potential which is 
the generalization of (21) 

_ * * 

U{ipi,ip2)^ li-lJ'livtvi)- tJ'liVl^'l)- IJ'1{V2¥>2)- fJ'l{<ptV2)] (37) 
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+A3((^]^(^l)((^^(^2)+ A3 {,(ptVl){<pt<P2) 
* 

+A6((^]^(^i)((^|(^2)+ Ae {^t^i){^2^i) 

* 

+A7((/7j(/72)(</?l"</'2)+ A7 {V2V2){vtV\) 

a denotes a complex conjugated. The potential terms with complex param- 
eters /x^2 ^5) ^6; ^7 explicitly violate CP invariance. With the help of 
U{1)y hypercharge symmetry of the model the phases 9^ and ^5,6,7 of com- 
plex parameters /i^2 -^5,6,7 can be removed by the phase rotation of scalar 
doublets <fi^2- In the case when vi in (3) is taken real and positive, and r) 
is the overall phase of ip2 doublet, the conditions to remove explicit phases 
from the potential are [20] 

61^ - 77 = n^TT, 65 -2r]^ n^TT, Oqj - rj ^ riejTr (38) 

where n^, rii {i =5,6,7) are arbitrary integer numbers. Or equivalently, in 
terms of complex parameters the conditions for the absence of explicit CP 
violation in the effective potential (37) are [5] 

* * 

lm{i^t^ As) = 0, Im{i^l^ Ag) = 0, lm{i^l^ A7) = 0. (39) 

The phase of complex /x^2 can be always rotated away, so /xf2 can be taken 
real. In the scenario of fine-tuning for the phases, when the conditions (38) 
or (39) are satisfied, Aj are also real. Otherwise Aj will be redefined after a 
phase rotation, keeping explicitly CP violating terms. 

If the phases of /i^2, \ a-re rotated away so there are no exphcitly CP 
noninvariant potential terms, CP invariance can be nevertheless broken spon- 
taneously. Using the convention when Vi and V2 in formula (3) are real and 
positive and selecting the ^ phase-dependent terms of (37), which arc of the 
form acos^,^ + 6cos^, where a = Asf ^^1/2 and b = {\Qvf/2 + X7V2/2 — nl2)viV2 
(0 < ^ < tt) one can find their minimum (see the appendix of [20]) at A5 > 
and 

,,2 1)2 Y 2 i>2 \ „2 

cose = " T ? ~ ^ ' (40) 
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(the special case yU^2 =0 "^^^ found in [1, 3]). If | /x^2 ~ T^e^^ ~ T^^"^! I — 
Asf^s^c^ then the spontaneously CP violating extremum does not exist and 
the minimum of the potential takes place at the endpoints cos^ = ±1. For 
6 > the minimum is at ^ = tt and for 6 < the minimum is at ^ = 0. 
If As <0, the extremum inside the interval |cos^| < 1 is the maximum, so 
the local minima are at the endpoints = 0, tt. The case ^ = is reduced 
to the case ,^ = tt by the change of sign for Ag,A7,/i^2 (change of sign for 
the b). The special case of CP conservation takes place at ^ = 7r/2 when 
/i^2 = T-^sC^ + Y^'^^p- ^^^^ purely imaginary < (p2 > [21] our di- 

agonalization procedure must be reconsidered. For instance, CP-even Higgs 
mass eigenstates are formed in this case by different orthogonal linear com- 
binations of real and imaginary parts of scalar doublets, Re(/9i and Imy92. 

The substitution of (26) to the no-extremum condition |/i^2 ~ Y-^^c^ ~ 
■yAysll > X^v'^spcp gives >0. In the case A5 <0 the absolute minimum of 

the potential takes place at ^ =0 (but not ^ = tt) if yU.^2 — T ^^c^ — ^ Ays^j > 0, 
which gives in combination with (26) m\ >\ A5 | v"^. Even if A5 is of the 
order of 1, the latter condition is valid when ttia is in the mass range of 
the order or greater than v. In the case of real parameters ^^2^ ^« ^'^^ 
straightforward to combine spontaneous CP violation with our procedure 
of diagonalization. In the case of complex parameters the situation may 
be changed. The extremum conditions can be found from the study of the 
fourth power equation with the coefficients depending on real and imaginary 
parts of /i^25 ^i- Nevertheless we are not going to consider spontaneous CP 
violation further on. With real parameters /x^2) ^« the absence of 

spontaneous CP violation the minimum of the potential (37) can be taken at 
^ =0. At the same time insofar as the physical motivation of the fine-tuning 
conditions (38), (39) for CP conservation is not available, in the following we 
consider the general case of diagonalization and minimization of the two- 
Higgs-doublet potential with complex parameters. 

For the diagonalization of (37) in the ground state we used the ansatz 
(22)-(28) to be taken for real parts of parameters. The real part of 
parameter is expressed through the real parts of X5,6,7 using (26) 

Re/x^2 = m^s^c^ ^^(s/jc^ReAs ReAe + ReAr) (41) 

defining also the real parts of Ai_2,3,4 and /if, by means of (22)-(25) and 
(27), (28). The substitution of complex //j and A^ to the potential (37) leads 
to the linear term and the non-diagonal mass term which are dependent on 
the imaginary parts of //i2) K 

U{ipi,(p2) = coA + cihA + C2HA (42) 
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Zi Zi Zi 

+ third and fourth order terms in /i, if, A, ii^ 

where 

ly^ — it' — — 

Co = -V Im/i^2 + yS/3C/3lmA5 + y (c^ImAe + s^ImAy) (43) 

Cl = ImAtL + ^(s2/3Sa-/3 - 2Ca+/3)lmA5 

—^{Sli(^p(^<x-fi - 3SaC^)lmA6 + —{SfiCfiC^^ji - 3c„S;3)lmA7 

C2 = Ca-/3 Im/ii2 + ^(c2/3Sa_/3 - 3s„+/3)lmA5 

— — 
-y + 2CaC/3)linA6 - y (s|ca-/3 + 2SaSp)lm\-j 

In the case of CP conservation (38), (39) the hnear and non-diagonal sec- 
ond order terms hA and HA do not appear because all imaginary parts of 
parameters can be removed. The linear term in A demonstrates that after 
the introduction of complex couplings we can be out of a local minimum of 
the potential U{<fi, ^2)- The minimization condition for the imaginary parts 
Co =0 must be imposed. If Aj, i =5,6,7 are zero (this occurs in THDM if 
we additionally introduce a global U{1)q symmetry, softly broken by the di- 
mension two /i^2 term [22]), the imaginary part of /xjg can be removed by the 
phase rotation of (/?2, so the tree-level potential is CP invariant. In the gen- 
eral THDM with nonzero parameters the phase rotation of (^2 which removes 
Im/i^2 redefines Aj, i =5,6,7. These simple observations for the THDM poten- 
tial isolated from any other physical fields are no longer trivial if, keeping in 
mind the MSSM, we switch on the interaction of ipi and <^2 with scalar quarks. 
CP invariance of the latter is (softly) broken by the dimension three terms 
with the Higgs mixing parameter /x and the trilinear parameters A^^i, of the 
form /i,2^6,tV'i 2'i*R^L, Ai/2,iV^i 2'i*L^R {1 = ii b; ipl 2 are the neutral components 
of the Higgs doublets, /i,2 = V^^fe,t/'Wi,2)- Then quartic scalar interaction 
parameters Aj, i =5,6,7 are affected by radiative corrections from one-loop 
diagrams with scalar quarks of the order of jp / Mgjj gy A'-^^/^lc/sy and 
jiA^ /MgjjgY [5], so the phases of A,, i =5,6,7 are defined by the phases of 
complex /i and A, thus constraining the phase of //^2 power of the condi- 
tions (41), (43). In the case of the Born level MSSM potential with a global 
U{1)q symmetry when Aj =0 (i =5,6,7) the complex /i^2 parameter can still 
appear beyond the tree-level due to the same CP violating Yukawa interac- 
tions of scalar quarks with the Higgs fields. This possibihty of the ii\2 phase 
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induced in higher orders by radiative corrections calculated diagrammatically 
has been considered in [23]. The restoration of potential minimum can be 
achieved by means of the opposite sign quantum correction term, originating 
from the tadpole diagrams with the pseudoscalar A connected to the squark 
loops ^. 

In the classical minimum cq =0 we find 

0,2 



Cl 



(44) 



C2 = -y(SaC| + CaS|)linA5 - t"^ (CaC/^ImAe + SaS/jImAy) 

The second order terms hA and HA in (42) can be removed as usual by the 
orthogonal rotation aij =1,2,3) in h,H,A sector 

/ h 
H 



(/i, H, A) M' 



(/ii,/i2,/i3) a^fc aij 




(45) 



where the mass matrix has the form 




C2 



Cl 
C2 



(46) 



Squared masses of the physical states hi, h2, h^, which are the Higgs bosons 
without definite CP-parity, are defined by the eigenvalues of mass matrix 
(roots of the cubic equation for eigenvalues are given by Cardano formulae) 







?)cos(-) - - 




(47) 


mil 




?)cos( g ) 


02 
"3 








e + An 
q)cos{ g ) 


02 
"3 





where 



r = — (9aia2 — 27ao "~ 203), 
22,22 222 



02 = 



-rriu 



m 



H 



arccos 



Q = g(3ai -al) 

22, 22, 22 2 2 
ai — myjrijj + rrifjn^ + mjjrrij^ — c-^ — Cg, 



^however, with nonzero A5, Ag and Ay, the factor of the scalar-pseudoscalar Higgs coun- 
terterm is not explicitly proportional to the tadpole renormalization constant, or the tad- 
pole parameter cq. 
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One can see that in the hmiting case of CP-conserving potential Ci,2 the 
following correspondence takes place: rUh, itt-h and m/jj — > 

TJiA- The normalized eigenvectors of the matrix M^, which are at the same 
time the matrix elements of a^, {h, H, A) = Qijhj, have the form a^j = a^j/rij, 
where 

= ii^H - ^hJi^A - ^IJ - 4), 4i = C1C2, a'^i = -ci{mjj - mlj 
' I 2 2 \ ' I 2 2 \ ' / 2 2 \/ 2 2 \ 



and rii — ^ (oii + 02i + '^3^) ^ • Representations for the triple and quartic Higgs 
boson self-interactions in the case of CP violating potential are given by the 
expansions of structures ttijhjaikhkauhi, and aijhjaikhkauhiai^mhm, they are 
bulky and not very telling, so we do not show them here. If the imaginary 
parts of Xq and A7 are not small, large off-diagonal elements of the mixing 
matrix Uij could appear leading to significant mass splittings of the Higgs 
states and modifications of the Higgs boson interactions. 

We assume that in the Yukawa sector < (fi > couples only to down 
fermions 

(^lYi — — 

Yl ^^^Trm ~ bPl{i+l5)'4^2a<fl+i^2a{'i--l5)Aft] (48) 

a=d,s,b 2V2mwSwCf} 

where {h, H, A) — Uijhj, for the first generation quarks V^i = {u, Vudd+VusS+ 
Vubb}, ip2a = {d, s, b) and analogous terms for c and t quarks {Vab denotes 
the CKM matrix elements), and < ip2 > couples only to up fermions (model 
of type n [24]): 

GTTl — — 

■[^1(1 + j5)ir2'ip2Vt + ^'2(1 - 75)«^2^lV^2] (49) 



where again physical hi,h2,hs states are introduced by means of the aij 
rotation, ipi — {u, Vudd + VusS -\- YuJ)}, ip2 — u and analogous terms for c 
and t quarks. 



^No ordering of masses to^i < 'mh2 < ™/i3 is required. If we want to keep this ordering, 
then aij written above, vahd for the case mn > m-A, must be changed. For the case 
niH < tha one should replace and change the sign of a^, Ois in the expressions 

for a,; 
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4 Higgs-gauge boson and Higgs-fermion cou- 
plings in the MSSM with expHcit CP vio- 
lation 

In the following we shall focus on the MSSM scenario for the two-Higgs- 
doublet model, which allows to restrict strongly the THDM parameter space. 
It is not the only one possible, Standard Model-hke scenarios in the general 
two-Higgs-doublet model have been discussed in [25]. Detailed consideration 
in the framework of MSSM has been performed in [5] (also [6]). In this sec- 
tion we would like only to compare qualitatively our results with the results 
of these approaches. Our calculation follows somewhat different scheme. In 
[5] the tree- level two-Higgs doublet potential is CP invariant. The phase ^ of 
//^2 is radiatively induced by the tadpole diagrams and can be absorbed in 
the definition of the /x parameter which appears in the stop mixing matrix 
off-diagonal element Af — fi/tgf3. The A5, Xa and Ay terms are also radiatively 
induced by the threshold effects. At the same time the trilinear couplings 
Af, Ah also carry a phase ^ , so both the radiatively induced and the trilinear 
phases contribute to the phase argdiA) of the A5, Ag and A7 terms. We do not 
account for the radiatively induced phase which is calculated diagrammat- 
ically. In the case under consideration the fine-tuning conditions (38), (39) 
are not fulfilled, so CP invariance of the two-doublet potential is explicitly 
broken by complex parameters at the tree level. The real and imaginary 
parts of //^2 parameter are defined by means of the condition (41) for the real 
parts of ^5,6,7 parameters and the minimization condition (43) (where Cq =0) 
for their imaginary parts. In the following calculations complex parameters 
A5,6,7 are specified in the framework of the MSSM. 

We used the two-loop symbohc results for Aj, i =1,...7, which were ob- 
tained in the RG approach [17] and extended to the case of CP violation in 
[5]. The parameters A5, Xq and A7 are nonzero in the next-to- leading order 
approximmation (RG improved leading order approximation), so using (41) 
and (43) 

- 1 - 1 
Re/i^2 = m\spCi3 + v'^{spCi3ReX5 -\ — ReAe H — ReAy) (50) 

v"^ - - - 

Im//^2 = -z-{si3CplmX5 + c'plmXG + s'plmXj) (51) 

where As^e,? depend on the finite term corrections to the leading logarithmic 
result which appear from the one-loop diagrams with trilinear couplings. The 

"^for a rewiev see e.g. [26] 
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analytical representation of [5] has the form 



1927r2 "MlusY 167r2 

|2 



96x2 "Mst/si- Ms[/sy M|[,sr" 16ir2^2 ' 2 



Atfi and /x are the factors in front of Higgs-squark(left)-squark(right) trilinear 
terms, Msusy is the SUSY energy scale, rritop-i mi, are the on-shell running 
masses of the third generation quarks, and t = log ^^^J^ , ht = ^^^""^ , 

hh = , gs = ^/4:7ras- The trilinear parameters At, and jj, of the Higgs 
boson interaction with the third generation squarks can be generally speaking 
complex. In this case the Aj, i —1,...7 parameters of the two-doublet Higgs 
potential are defined by tg/9, SUSY scale Msusy, and six relevant parameters 
in the sector of Higgs boson interaction with the third generation squarks: fi, 
arg(yu), At, arg(At), Ab, arg(A5). In the following consideration for simplicity 
we assume \At\ = \Ab\ and assign the universal phase 9 to iiAt and nA^ so 
that 9 —a,-rg{ij,At) =aTg{fxAb). Then using the exphcit structure of (52) CP 

invariance conditions (39) can be rewritten in the form Iiii(a*^2 A*^) =0 [5]- 
The couphngs of W and Z bosons to the hi, /i2, hs scalars have the form 

Vij. V„ hi fv9p.u{ca-pa2i - Sa-paii) 

Vfj, K h2 fvgixACa-130'22 " Sa-fSa^) 

Vn K ^3 fvgiiu{Ca-l3a2Z " s^-paiz) 

where V ^ W,Z, fw = ±:mw and fz = —^mw The couplings of 
hi,h2, hs bosons to the t and b quarks have the form 
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ftj-{Sa0'21 + Caflll - ^0^03175) 
/< — (S„a22 + Caai2 - iCpa32l5) 
ftj^{Saa23 + Caai3 - iCfsass'J^) 

fb — {Caa23 - SaOlS " ^5/303375) 



The Higgs boson mass spectrum of the CP conserving hmit 6 =0 (in 
this hmit ttij = diag{l, 1, 1}) is shown in Fig.l. For the case of exphcit CP 
violation in the two-doublet Higgs potential we take the parameter set pL —-2 
TeV, At = Ah =-1.8 TeV, Msusy =0.5 TeV, =220 GeV, tg/3 =4 which is 
typical for the region of MSSM parameter space where the imaginary parts of 
A5, Xq and A7 are large (of the order of 0.1-1) ^ . We demonstrate in Fig. 2 the 
neutral Higgs boson masses given by (47) and the mixing matrix elements Uij 
as a function of the universal phase 9 —ATgdiAt^i,) . The Higgs boson masses 
of the CP conserving limit are substantially changed when the phase 6 is not 
small. The rrihi in Fig. 2 is always smaller than ruh and ruhz a downfall at 
the phase values around 7r/4. The Higgs- vector boson WWhi, ZZhi and the 
Higgs-fermion qqhi (g=t,b) interaction vertices as a function of the phase 9 
are shown in Fig. 3. One can observe that the hi couphngs to gauge bosons 
W, Z decrease by about 15% if the phase of A5, Ag, Ay is large enough. 
Nonzero couplings of to gauge bosons appear. The changes of the bbhi 
and the 66/12 coupling regime are also rather pronounced (see Fig. 3). In the 
region of MSSM parameter space where the ruhs is around 150-250 GeV and 
the II and At^i, parameters are of the order of TeV the regime of strong mixing 
in the Higgs sector takes place. As a result the light Higgs boson hi could 
have not been observed at LEP2 {^/s =200 GeV) in the production channels 
e+e^ — > hiZ, e+e^ i^e^ehi for the reason of suppressed couplings to gauge 
bosons, while the /i2,/i3 bosons are sufficiently heavy to be not produced 
on mass-shell at the LEP2 energy. Detailed analysis of this scenario can be 
found in [5. 6]. 

detailed discussion of possible combined constraints on the MSSM parameter space 
from the cosmology, direct searches and indirect measurements (rare decays) can be found 
in [27] 



where ft^ 



i t hi 
i t h2 
1 1 /is 

hhhi 
6 6/i2 
6 6/i3 

e 

2si4/ mw 
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5 Triple and quartic Higgs boson couplings 
in the MSSM with expHcit CP violation 

In the regions of the MSSM parameter space where the couphngs of hghtest 
Higgs boson hi to gauge bosons and top quark are suppressed, traditional 
channels of Higgs boson production by radiation from W, Z or t line and 
WW, ZZ fusion can have too small rate to be experimentally observed. For 
this reason it is interesting to consider the possibility of double Higgs pro- 
duction (like gg ^ h2 ^ hihi) defined by the self-coupling vertices. Such 
calculations are known in the CP conserving limit [8] , when the cross sections 
of double and triple Higgs boson production turn out to be very small. Only 
some of them are accessible for observation at high luminosity colliders. In 
the case of CP violation some self- couplings can be substantially increased, 
providing better possibilities for the experimental reconstruction. 

The A5, Xq and A7 potential terms can modify significantly the Higgs 
boson self-interaction vertices calculated in the leading one-loop approxima- 
tion with Aj =0 (i =5,6,7). At the next-to-leading order approximation the 
Aj couplings (52) include the terms of the order of h^^^'^A^^^/M^jjgy and 
h^^^liAi^i,/ MgjjgYi so they can reach the values of the order of 0.1-1 at moder- 
ate values of Msusy and jd, and At^b taken at TeV energy scale . For example, 
in the CP conserving limit 6 =0 the hhh vertex in the mass parametrization 
has the form 

3e 

ghhh = [-{cpcl- spsDml + cl_^cp+am\ (53) 

mwsw'^20 

The contributions of As^g,? terms and the m\ ^ terms in this expression are 
of the same order if As^gj ~ The rotation of the 9 =0 mass eigenstates 

by the matrix defined by (45) gives the ghihihi vertex of different form 
but also with substantial contributions of the A5,6,7 terms. 

Using the parameter set described in the previous section, we show the 
values of various triple and quartic Higgs boson self-interaction vertices as a 
function of the universal phase 9 =axg{^At) =aig{iiAi)) in Fig. 4. The values 
of Higgs boson self-interaction vertices in the CP-conserving limit 9 =0,7r 
and in the leading order approximation A5 = Ag = A7 =0 are marked in Fig. 4 
by horisontal arrows. The A5, Xq and A7 potential terms induced in the next- 
to-leading order approximation introduce very large corrections to the triple 
and quartic self-interactions of Higgs bosons. In the region of the MSSM 
parameter space under consideration the difference of the leading order and 
the next-to-leading order vertex factors can be several times in some ranges 
of the phase variation. 
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6 Summary 



We demonstrated the tree-level equivalence of the two-Higgs-doublet model 
potentials (1) and (9), where CP-invariance can be explicitly broken by the 
Xq term in (1) or by the complex /i^2) ^5 terms in (9). The parameters Aj 
(i=l,...6) of (1) and fif, ^^2, Aj (i=l,...5) of (9) arc related by the equations 
(10). In the case of real parameters the diagonalization of potential (9) in 
the ground state can be performed by means of the substitutions (14)- (20) 
which express the A, and //^, /i| parameters through the Higgs boson masses 
"iT^h, TTT'H, TTT-A, '"^i/ij thc mixing angles a, [3 and the /i^2 parameter. In the 
general case the Ag and A7 potential terms (21) should be also considered with 
the diagonalization and minimization conditions (22)- (28). If the complex 
parameters /x^, /Xg; (i=lr--7) and /if 2 are introduced, CP invariance of the 
hermitian potential (37) is exphcitly violated at the tree-level unless the fine- 
tuning conditions (38) or (39) for parameters arc satisfied. So in the following 
we consider the problem of diagonalization in the local minimum for the two- 
Higgs-doublet potential which is not CP invariant. For the diagonalization 
of potential (37) again we use the substitution (22)-(28) to be taken for real 
parts of parameters. The minimization of potential (37) at the tree level takes 
place with the condition cq =0 (43) for the imaginary parts of parameters. 
The imaginary parts of A5, Ag and A7 give rise to the CP-odd/CP-even Higgs 
boson off-diagonal terms, which are removed by the orthogonal rotation in 
(/i, i/. A) space, giving mass eigenstates hi,h2, without definite CP-parity 
and with the mass spectrum and couplings substantially different from the 
masses and couplings of CP-cvcn and CP-odd states h, H, A, if the imaginary 
parts of parameters A5, Ag and Ay are sufficiently large. 

In the framework of MSSM the real parts of Aj (i=l,...5) couplings are 
fixed at the SUSY energy scale by the conditions (29). Radiative corrections 
to the Xf^^^ (i=l,...7) couphngs are generated at the mw energy scale. The 
equations (31)- (33) express the mixing angle a and masses of Higgs bosons in 
terms of the radiative corrections to Af^'^^ (i=l,...7) couplings (e.g. given by 
the RG evolution). They are valid independently on the particular scheme 
which is used for calculation of radiative corrections to the Xf^^^ (i=l,...7). 

In the next-to-leading order approximation the complex A5, Ag and Ay 
parameters are generated by the soft CP violating Yukawa interactions of 
Higgs bosons with the scalar quarks. Using the results of [5] we calculated the 
Higgs-gauge boson, Higgs-fermion and the Higgs triple and quartic couplings 
for a representative MSSM parameter set, when the off-diagonal elements of 
the Higgs boson mixing matrix are large. The A5, Ag and Ay parameters in- 
troduce significant corrections to the Higgs self-interaction, even in the case 
when their effects on the Higgs-gauge boson and Higgs-fermion couplings 
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are rather small. These corrections could rather strongly (by one-two orders 
of magnitude in comparison with the case of CP conservation) enhance or 
suppress some channels of multiple Higgs boson production at next collid- 
ers, providing discriminative tests of CP violation in the Higgs sector and 
improved feasibility to reconstruct experimentally the Higgs potential. 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


h h h 
H H H 
H H h 
H h h 
H A A 
h A A 
h H+ H- 
H H+ H- 


M^Zsd ^2/3(c^C/3 sis^)mi + 2ci_^Co,+^^ii2] 

L sl^ [s2/3(Sa4 + Ca4)"^H + 4pCa-fim\ 2s«+;3/i?2] 

2p 



Table 1. Triple Higgs boson interaction vertices in the general 
two-Higgs-doublet model, /X12 parametrisation. 



Fields ill the \'ertex 



Variational (leii\'ati\'e of Lagrangian b}' fields 



h h h 

H H H 

H H h 

H h h 

H A A 

h A A 

h H+ H- 

H H+ H- 



Mw S23 
3e I 

2mJ::'.,, [-(2^1/ + ^Ds2a + 2(3s,Cc, + S/3C^)(mi + i;2A5)] 



2Mw Sw S2fj I 



e 



Table 2. Triple Higgs boson interaction vertices in the general 
two-Higgs-doublet model, A5 parametrisation. 
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Fields in the vertex 


Variational derivative of Lagrangian by fields 


h h h h 
H H H H 
AO AO AO AO 
H H H h 
H h h h 
H H h h 

H H A"^ 

/i /i AO AO 

AO AO /i 

//+ H+ H- H- 
H+ H- AO AO 
H+ H- h h 

H+ H- H H 

H H+ H- h 


2 

4 M^' [^S2l3{clcp Slsp)ml + S2l3S2aSa-pm]j ^Ca+l3\lA2 

4 s2 „3 [ S2/3S2a(3S2aSQ,_3 4:S a-^ pCa^ i3 2Sa+i3Sa-l3Ca-j3)'n^h 

W 2/3 

+S2/3S2a(s2/3 + SS2aSl_p)ml - 8{3slcl - 44)/^?2] 
4M^4s^^[ '^^2pS2aSa-l3(CaCp So,S^)m^ 2s2^C„_^m^ 
-S2/3(s2aS2/3 + 3s2_^s2^^ _ S2/3S^_^)m|^ + 4:{4pSl_p + S^+^)/X?2] 
4 M^fei; [ «i^(4c2/3C2a + 3s^_^s2_^^ ^ S^_^)m| 2si^s2 _^^2 

-2S2/3S2aCa_^(s«4 + Ca^)"^!^ + 2(si/3Sa-/3 + 4(Ca4 ~ ^<^^W)l^l2\ 
4M^s|,s^^[ 2S2/3S2aCa-/3(CaC| Sa^)"^^ + *2/3Sa-/3Ca-/3"^^ 

-2S2/3S2aSc,_/3(Sa4 + (^aS^p)^H + 2(2S2aC2^ " S2/3S«_^C^_/3)/U?2] 
~^m2,4 s^^ [g2/?(Ca4 - SaS^p)n^h + S2p{SaCp + Cagp)m]j - 2ci^//f2] 

4 M2,4.3^ [ S2^(4c2aC2;3 + ^sl_psl^p + S^_^)m| 2slpsl_f^m]j± 
-2S2pS2aCa-l3{SaCp + CaS^^)mjj + 2(s|^S^_^ + 4(c«4 " Sa^) V12] 
4M^s^5^^[ 2S2/3S2aSa-/3(Ca4 ^a^)"^^ 2s^^4_^m2^± 
+S2/3(s2aS2^ - 3s2_^s2^^ + si_p)m\ + 4(cl^S^_^ + s2^^)/X?2] 
2 si [ ^2pS2aCa-p{CaCl3 SaS p)mf^ + S2i3Sa-pCa-pmjj± 
-S2f3S2aSa-p{CaS$ + SaCl)mjf + 2{2s2aC2f3 - slffSa-pCa-fj) HI2] 



Table 3. Quartic Higgs boson interaction vertices in the general 
two-Higgs-doublet model, /ii2 parametrisation. 
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Figure 1: Masses of the neutral and charged Higgs bosons h,H,H^ versus 
the pseudoscalar mass and the trihnear constants At, calculated by 
means of (31)-(33) with the analytical Aj (i=l,...7) parametrization of [5]. 
The AA5 is chosen to be positive. The CP-conserving limit 6 =0 is taken, 
(a) tg/3 =4, MsusY =0.5 TeV, At ^ A, ^ pt =0; (b) tgP =4, Msusy =0.5 
TeV, At = Ab =0.9 TeV, /j, = -1.5 TeV; (c) tg/3 =4, Msusy =0.5 TeV, 
niA =220 GeV, =0, At = At; (d) tgp =4, Msusy =0.5 TeV, tua =220 
GeV, fi = —2 TeV, At = Ah. Very small variations of the charged Higgs 
boson mass mH± in (d) are due to the cancellation of leading power terms 
~ IJ^Al^/MgjjgY, see [5], in the difference of AA4 and AA5, see (33). If 
AA5 is chosen to be negative (purely imaginary ^A in (52)), mu increases in 
comparison with the case At — A^, — ^ =0. 
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Figure 2: Masses of the neutral Higgs bosons and the mixing matrix elements 
as a function of the Xq and Ay phase. The Aj parameters are taken from [5] 
at the MSSM parameter values tg/3 =4, tua =220 GeV, Msusy =0.5 TeV, 
At = Ab = -1.8 TeV, = -2 TeV. 
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Figure 3: Higgs-gauge boson and Higgs-fermion vertex factors as a func- 
tion of the Ae and Ay phase. The Aj parameters are taken from [5] at 
the MSSM parameter values tg/3 =4, tua =220 GeV, Msusy =0.5 TeV, 
At — A}, — —1.8 TeV, ji — —2 TeV. For the couphng with fermions we plot 
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Figure 4: Triple and quartic Higgs boson vertex factors as a function of the Ae 
and A7 phase. The Aj parameters are taken from [5] at the MSSM parameter 
values tg/3 =4, tua =220 GeV, Msusy =0.5 TeV, ^ = A = -1.8 TeV, 
/i = —2 TeV. Horisontal arrows indicate the values of vertex factors in the 
CP-conserving limit 9 —0 and the leading order approximation A5 = Ae = 
A7 =0. 
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